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We investigate the properties of cuprate superconductors subject to apphed current, using modi- 
fied Gutzwiller projected d-wave BCS states. The parent states include quasiparticle and quasihole 
pockets, of variationally determined size, generated by the current. We identify two different mecha- 
nisms for the destruction of superconductivity at the critical current: at high hole doping (x ^ 0.15) 
the pockets grow and completely destroy the gap, in a BCS-like mechanism; in the underdoped 
regime, the superfiuid stiffness vanishes at a maximal phase twist with pairing still intact. This 
result is indicative of a pseudogapped 'normal' state which retains pairing correlations. The critical 
current as a function of doping displays a dome shape, similar to Tc. We predict unique signatures 
of the current induced Fermi pockets that can be seen in angle resolved photo emission spectroscopy. 

PACS numbers: 74.72.-h, 74.25.Sv, 74.20.Mn 



Introduction - Shortly after the discovery of high 
temperature superconductivity in the cuprates, Ander- 
son proposed that the strong on-site repulsion and prox- 
imity to the Mott insulator play a crucial role in deter- 
mining the electronic correlations in these materials [l[. 
Such correlations are naturally included in Gutzwiller 
projected variational BCS wave- functions, which indeed 
provide quite a good description of the superconducting 
state [1, 3]- However, variational studies have not ad- 
dressed the key open questions concerning the destruc- 
tion of superconductivity with temperature, the nature of 
the seemingly non Fermi- liquid normal state, the pseu- 
dogap phenomenon, possible competing orders, and an 
associated quantum critical point. 

Here we take a step toward addressing these issues us- 
ing Gutzwiller projected states. The idea is, rather than 
investigate the effects of temperature, which are inacces- 
sible in this approach, to study the consequences of an 
applied current at zero temperature. Both mechanisms 
similarly lead to the eventual demise of superconductivity 
and the establishment of a normal state. The reduction 
of the superfiuid stiffness ps{T) is linear in temperature 
at low T due to thermally excited quasiparticles in the 
(i-wave gap nodes [4]. A similar linear reduction is ex- 
pected to occur with external-current due to formation of 
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FIG. 1: The low energy quasiparticle spectrum of a d-wave su- 
perconductor with a Doppler shift due to a superflow of wave- 
vector Q = Qx. The ground state consists of the marked 
quasiparticle and quasihole pockets. 



quasiparticle (QP) pockets around the nodes, which are 
Doppler-shifted away from zero energy in the presence of 
superflow (see illustration in Fig. [1]). Our approach 
can therefore lend insight into the puzzle concerning the 
evolution of the slope dps/dT with hole doping [gI, 0- 
Furthermore, the behavior of the superconductor with 
uniform current is intimately connected to the effects of 
a magnetic field, which generates supercurrents around 
the vortex cores. The critical current which we calculate 
is related to the maximal circulating current found at the 
edge of a vortex core and the state formed at the critical 
current is connected to the normal state in the core. 

Whether superconductivity is suppressed by temper- 
ature, uniform current or (orbital) magnetic field, there 
appears to be a qualitative change in the mechanism of 
its destruction and in the corresponding 'normal' state, 
upon doping the system from the underdoped to the over- 
doped regime A common view is that in the under- 
doped regime superconductivity is destroyed with pairing 
still intact [8, 9], giving rise to a pseudogap above Tc [Toj . 
In the over doped regime by contrast, the transition is 
perceived to be conventional BCS-like. We shall directly 
test such hypotheses using the variational approach by 
suppressing superconductivity with a current. 

To allow for the presence of a supercurrent we extend 
the commonly used family of Gutzwiller projected d-wave 
BCS variational states. The current is embodied in Fermi 
pockets containing the quasiparticles and quasiholes gen- 
erated by the Doppler shift around the gap nodes. A new 
variational parameter determines the size of the pockets 
for a given superflow. We use these variational states 
to compute the critical current and discern the different 
mechanisms that lead to destruction of superconductivity 
at varying hole doping, as described in the abstract. 

The variational wave-function - We set out to 

investigate the ground state of the t-t'-J Hamiltonian on 
a square lattice subject to superflow. In order to later 
fix the superflow wave- vector Q = Qx it is convenient to 
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FIG. 2: Discontinuities in (rik) a?^c? weight, (a-b) Gray- 
scale plots of Zk = I Vk^kl on a 24x24 lattice, doping x = 0.12 
and Q = (a) and Q = 7r/2 (b). (c) The QP weight on the 
current-induced Fermi surfaces in a d-wave BCS state, (d) 
QP weight Zk((5 = ^/2) (black circles) and the QP weight 
renormalization ak = Z\^/Z^^^^ (squares) along the quasi- 
hole pocket Fermi surface (see inset), (e) Fermi surface av- 
erage of ak as a function of doping (squares) compared with 
the zero current result of Ref. [3] (dotted line) and the renor- 
malized mean field result = 2x/(l + x) (dashed line). 



impose a constant vector potential along the x direction 
via a Peierls substitution, 



Here Pq = 11^(1 — n^^n^^) is the Gutzwiller projection 
operator. The variational wave- functions we use are 
Gutzwiller projected d-wave BCS states containing filled 
QP and quasihole (QH) pockets around the nodes. Be- 
sides the usual parameters A and /i we include a new 
variational parameter Qvar which determines the size of 
the current-induced pockets. The superflow is imposed 
by constraining A to be spatially uniform (and real) re- 
gardless of Q. 

To motivate the construction of the wave-function, it 
is useful to recall the mean field BCS description of a 
d-wdive superconductor subject to current of wave- vector 
Q. The Bogohubov De Gennes spectrum is given by 



E^iQ) = ?-,k(Q)±^Ci,k(Q) + A^ 

« ^VF(k)-Q±^k(0) 



where is the single particle dispersion, Ak = 
A(Q)(cos kx—co^ ky) the d-wave gap function, and (f±^k = 
(^k+Q/2 i ^k-Q/2)/2. The second line above gives the 
low energy spectrum for small Q. The BCS ground state 
is defined by occupation of all the negative energy QP 
states. Because of the Doppler shift vp •Q/2, this implies 
filling an elliptical pocket of QPs around the two down- 
stream nodes and QHs around the upstream nodes, cf. 
Fig. [TJ More generally, due to interactions, the Doppler 
shift may be renormalized to avp • Q/2 = vp • Qvar/2. 
The parameter a is commonlyreferred to as the effective 
QP current renormalization [3, In the parent BCS 
state we shall allow the freedom to determine the size of 
the pocket via the additional variational parameter Qvar 
independent of the externally applied twist Q. 

The variational state we construct is a Gutzwiller pro- 
jection of the mean field ground state with N electrons, 

I [A, /i, Qvar]) =PgPn n 7L|^BCS(Q)). 

Cr,kG'P(Qvar) 

P/v projects the state onto the N-particle subspace, 
7^ (Qvar) denotes the wave- vectors inside the QP pock- 
ets, and I ^BCs(Q) ) is the vacuum of the QP operators. 



(2) 



^2)1/2 



With ul{Q) = l-vl{Q) = - 

We minimize the energy as a function of the variational 
parameters A, ji and Qvar using variational Monte Carlo 
(VMC) 0, with J = t/?> and t' = -t/4. The opti- 
mized energy and wave-function are then used to study 
the effects of strong correlations on a current-carrying 
d-wave superconductor. 

Quasiparticle properties of the projected state - 
We turn to investigate basic properties of the projected 
current-carrying state as would be observed in angle 
resolved photoemission spectroscopy (ARPES). Specifi- 
cally we consider (nk), related to the spectral function 
via the sum rule (n^) = Ea(4a^ka) = J^jiujAQ{k,uj). 
From the variational state we can extract discontinuities 
in (nk) which correspond to the QP weights at uj = 0. 

To study the effects of strong correlations we compare 
the structure of (nk) in the projected state with that 
of the underlying BCS wave- function. Within BCS, the 
T = spectral function has the form 



^BCS 



(k,c^) 



S[u 



Km 



(1) 



Consequently, (nk) has discontinuities of magnitude u'^ 
(vl.) along the QP (QH) pocket Fermi surfaces, as de- 
picted in Fig. [2fc). At Q = the pockets shrink to 4 
nodal points with a QP weight of 1 along the nodal di- 
rection at each point [12|. 

For the comparison with BCS, we plot |Vk^k| = 
in the projected state with Q = and Q = 7r/2 x 
[Fig. [2ja,b)]. The pocket Fermi surfaces in Fig. EJb), 
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FIG. 3: Critical flow vs doping calculated from VMC (solid 
lines) vs RMFT (dashed lines) results, (a) Critical flow wave- 
vector Qc- (b) Length of residual paired region in the under- 
lying Fermi surface, (c) Critical current Jc vs doping. Right 
panels: pockets at the critical flow for three doping levels 
marked in panel (b). The pockets grow with doping until 
they touch the zone boundary at x > 0.15. 



marked by the peaks of Zk, behave similarly to the BCS 
case shown in Fig. Efc). To quantify the relation we ex- 
tract the ratio between the QP weight after projection 
and the bare BCS value along the hole-pocket boundary 
<^k = ^k/^k^^ [filled squares in Fig.[2fd)]. Interestingly, 
ak is almost independent of k. This is consistent with 
the central assumption of renormalized mean field theo- 
ries (RMFT) [13, Q, that the Gutzwiller projection can 
be encapsulated in a global renormalization factor. We 
also find that the renormalization factor is independent 
of current, and its doping dependence (Fig.[2fey) matches 
with the zero current result obtained in Ref. [3[. 

The above predictions for ARPES do not depend on 
the pocket shape. For clarity, we chose A, fi and Qvar 
to obtain circular pockets rather than the variationally 
determined elliptical shape. 

Critical current in the variational state - The 
critical current is the maximal current that can be im- 
posed on a superconductor 



dE 

Jc = max — - 
Q dQ 



(3) 



and it occurs at a critical flow wave- vector Q = Qc The 
maximum is taken under the constraint that at the flow Q 
the variational Fermi pockets do not yet cover the under- 
lying large Fermi surface, at which point superconductiv- 
ity would be lost due to complete depairing. Technically, 
we obtain Qc and Jc by fitting the variational energy, 
calculated on a 18 x 18 lattice, to a polynomial [l5| 



E{Q) = Eo 



0{Q% (4) 



The results of this analysis are shown in Fig. [3l Fig- 
ure [31(a) displays the variationally determined Qc as a 
function of the doping x. Notably, Qc is essentially in- 
dependent of X for X < xo = 0.15, and it drops sharply 
for higher doping. A complementary viewpoint is given in 
Fig.[3jb) showing q^^.^, the length of the residual gapped 
region of the Fermi surface, at the critical current. As 
doping is increased, the maximal pocket size grows and 
q . drops, until it reaches zero at x = xq. At x > xq 

• pair 

the critical current is obtained when the wave-function is 
fully depaired. The Fermi pockets at the critical current 
are shown in panels I-III in Fig. [3] for three doping levels. 

We conclude that optimal doping marks a transition 
between two distinct mechanisms for the destruction of 
superconductivity with current: A BCS-like mechanism 
for x>xo, and a "bosonic" mechanism, whereby super- 
fluidity is destroyed with pairing still present, for x<xo. 
The doping independent value of Qc suggests an analogy 
with lattice bosons for which superfluidity is destroyed at 
Qc = 7r/2 [l^. The residual pairing at the critical current 
suggests that the normal state formed upon destruction 
of superconductivity in the underdoped regime retains 
pairing correlations. 

The transition between the two mechanisms gives rise 
to a maximum of Jc at x a;o[l7|. The dome shape 
of the critical current as a function of doping, shown in 
Fig. EJc), is reminiscent of the observed doping depen- 
dence of the critical temperature, with a similar optimal 
doping level. We can estimate the expected critical cur- 
rent in the underdoped regime as Jc ~ psoQc- Taking 
PsO = pah from experiments, such as ^ and Qc ^ 7r/2 
gives a (3D) current density of Jc/d ^ lO^A/cm^, where 
d is the interlayer spacing. This is an estimate for the mi- 
croscopic critical current, or the maximal current at the 
edge of a vortex core. The macroscopic critical current 
may be significantly suppressed due to inhomogeneity. 

Before proceeding we note, that similar behavior might 
have been expected within RMFT [l3| , due to the compe- 
tition between the rising stiffness and falling gap which 
was noted early on [3]. Nevertheless, a calculation of 
the critical flow within RMFT (dashed lines in Fig. [3j) 
draws a markedly different picture than does the varia- 
tional calculation. In particular the critical flow wave- 
vector drops continuously as 1/x^ and the pocket size is 
always smaller than the Brillouin zone, so that there is 
no complete depairing at any doping level within RMFT. 

The current dependent superfluid stiffness - 

The current dependent superfluid stiffness is the non- 
linear response ps{Qo) = ^q^|q=Qo calculated 
within the variational scheme by fitting the optimized 
energy to Eq. ([4]) in the range 0<Q<Qc [lH. 

The result is displayed in FiglH^a), showing the dop- 
ing dependence of the variational stiffness pso (solid cir- 
cles) compared with the mean field calculation (dashed). 
Note that we obtain an improved estimate of the su- 



4 



0.05 

0.025 
0.075 



.(a) 




















Pdia 






^0.06 0.08 0.1 0.12 0.14 

doping (x) 



FIG. 4: (a) Zero current superfluid stiffness pso vs hole dop- 
ing in the underdoped regime (circles), compared with the 
diamagnetic response (squares), computed without allowing 
for pockets, and with the RMFT result for pso (dashed), (b) 
The slope 77 = dps/dQ vs doping. Compare with the much 
weaker doping dependence in the RMFT result (dashed). 



depends only weakly on doping, reflecting the puzzlingly 
weak doping dependence of dps/dT in experiments. 

Though our results pertain to the destruction of su- 
perconductivity with current at zero temperature, they 
are strongly suggestive of a pairing mechanism under- 
lying the pseudogap phenomenology in the underdoped 
cuprates at T > Tc. Further insight into the nature of 
the normal state can be gained by extending our work to 
allow competing orders, such as spin or charge-density 
wave, possibly reducing the variational energy as super- 
conductivity is suppressed 20|, Hi | . This may also eluci- 
date the role played by a quantum critical point, where 
the competing order in the normal state vanishes [22 1. 
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Foundation and the Minerva foundation. 



perfluid stiffness, including paramagnetic contributions, 
compared with the diamagnetic upper bound (squares) 
that was calculated previously [l8| . 

More importantly, we can now study the suppression 
of superfluidity with the flow, which at small Q has the 
form ps{Q) — psO — v\Q\- Figure lU^b) displays the coefli- 
cient 77 as a function of doping extracted from the varia- 
tional calculation (solid circles). This doping dependence 
is clearlyweaker than the RMFT result (dashed) in which 
T] ^ [5| . We note that in the underdoped regime 77 can 
be related to Qc via (j4]) as 7^ :^ pso/Qc- Since pso ~ x 
, the weak doping dependence of r] can be traced back 
to the doping independent critical flow Qc Hence it also 
stems from the bosonic mechanism of destruction of su- 
perconductivity in the underdoped regime. 

It is plausible that the weak doping dependence of 
dps/dQ shares the same origin as the unexplained weak 
doping dependence of dps/dT j6|. The latter is not un- 
derstood within mean fleld approaches 0, nor from cal- 
culations of the current carried by a Gutzwiller projected 
single QP state [ll^ . Both slopes stem from paramagnetic 
current of QPs near the nodes, induced either by temper- 
ature or current. Measurements of the stiffness reduction 
in the presence of supercurrent may clarify the relation 
between dps/dT and dps/dQ in cuprates. 

Summary - Using a new class of variational states 
we identified a transition between two mechanisms for 
the destruction of superconductivity with current. In 
the overdoped regime {x > 0.15) we found a BCS-like 
mechanism, where the critical current is reached due to 
complete depairing, while in the underdoped regime we 
observed a bosonic mechanism where superconductivity 
is destabilized while pairing survives. The critical current 
has a dome shape with an optimal doping set by this 
transition [l9|. We computed the linear reduction of the 
superfluid stiffness at small current. The slope dps/dQ 
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